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In a previous paper * the author has discussed the effect of viscosity in 
dissipating the energy of non-turbulent tidal flow of a sea, and the present 
one contains a simple extension of the method to non-periodic motion. In 
practically all parts of the open ocean there exist currents whose directions 
are almost constant and whose velocities may therefore be supposed to 
contain a component which is independent of the time. The motion is 
usually slow, the current seldom exceeding one or two knots, and on the 
average much less, so that in all probability the effect of turbulence will not 
be large. The work is therefore carried out on the assumption of a smooth 
flow, and a first approximation only is made. The results are used to obtain 
an estimate of the mean rate at which energy is dissipated in all the oceans 
by these permanent current streams. 

The second portion of the paper is devoted to a simple method of deducing 
the magnitude of the surface current solely from a knowledge of the 
distribution of the salinity or of the temperature of the water. There are 
not sufficient data available to obtain exact figures, but rough estimates of 
the orders of magnitude lead to concordant results. 

The Dissipation of Energy in Permanent Ocean Currents, 

Consider a portion of the sea so small that in the state of relative 
equilibrium its surface may be regarded as flat, and let it rotate about the 
vertical with an angular velocity o> which is to be taken as the component of 
the earth's diurnal rotation. Let the axis of rotation be taken as the axis 
of z, the origin being at the bottom of the sea and the equation of the mean 
surface being z = h. Take axes of x, y, rotating in their plane with the 
prescribed angular velocity &>, and let U, V, W be the components referred 
to these axes of the small relative velocity of the particle of water which at 
time t is at the point (x, y, z). 

* ' Roy. Soc. Proc.,' A, vol. 93, p. 348 (1917). 
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For the horizontal motion the equations are* 

at ay 



(1) 



(2) 



where the density is assumed uniform, v is the kinematic coefficient of 
viscosity, A is the Laplacian operator, and Z — H/g is the elevation of the 
surface above its position of relative equilibrium, II being the disturbing 
potential. 

Let U = U + lfUr cos (at — a) 

V= Y+Xv a cos (<rt-/3) 

where U, V, v, 9y v 9 are independent of t, and the summation includes terms of 
different periods. The parts of the solutions which are periodic functions 
have been discussed previously! on the assumption a and co are of the same 
order of magnitude and that a — 2 co is positive. By a very slight modifica- 
tion the work can be extended to the cases in which a is less than 2 co, that 
is, to the longer period motion. We proceed to consider the parts U, V of 
the solutions which are independent of the time, assuming that co is not zero. 
By equations (1) and (2) these are given by 

~2a>V = 



dx ' cb 2 



2*)U 



az , a 2 v 



s 



(3) 



where Z is the part of Z independent of the time, and A (U, V) have been 
replaced by 3 2 (U, Y)/dz 2 . 

Combining these equations we obtain 



a 2 -^ 



§T^(U±<V) 
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J v v \d% dy 

of which the solutions, subject to the conditions U, V zero at the bottom 
z = 0, and dTJJdz, dY/dz zero at the surface z = h, are 

cosh {7 (h—z) (1 ± l)} 



Zoo _ 



/_3 , _3 \ ^ 

\dx~ dyj J? 



(4) 



cosh {7^(1±0} 
where 7 2 = 00/ v. 

In the open sea yh is always very large, so that we can replace these 

solutions by 

U±*V = (U ±*Vo) [l-e-y*^] (5) 

where U 0) Vo are the surface values of U, Y. 

# See, for example, Lamb, ' Hydrodynamics,' §§ 2C6, 316. 
t * Eoy. Soe. Proc./ A, loc. cit. 
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Now let us suppose that Q is the local magnitude of the non-periodic part 
of the surface current relative to the earth's rotation, and that the local axes 
are chosen so that Q is in the direction Ox, Then Uo = Q, V = 0, and 
equations (5) become 

U = Q(l— 0~~Y*COS72) 

V = Q6~Y"sin7£ 



It is now easy to justify the neglect of the terms in obtaining equations (3). 
The non-periodic terms in the equation of continuity lead to the relation 

which from equations (4) is equivalent to 



f 



K-+3- )dz = 0, 
o\cte ay I 



Hence from the same equations 

liH?) (u ' v) = ' 

so that the neglected terms are, to this approximation, actually zero. 

The total rate of dissipation of energy per unit surface area of the sea may 
be taken* as 



->j> 



* / 3 2 U —B 8 ^, 



where p is the coefficient of viscosity. 

Substituting from equations (2) and taking the mean value over a long 
period of time we obtain a result which can be expressed in the form 



^fo( U "3? +V ^') &+2 ^ 



where E^ is the contribution from the terms of period 27r/cr. The mean rate 
of energy dissipation due to each term in the current velocity can thus be 
evaluated separately and the several results added. The contribution from 
the non-periodic terms given by equations (6) is found to be 7/7Q 2 to the 
order of approximation already used, thus depending only on the velocity and 
not involving its gradients. 

Writing co = X2 sin X } where X2 is the earth's diurnal rotation and X is the 
latitude of the place, we obtain for the mean rate, F, of dissipation of energy 
per unit surface area due to the permanent ocean currents 

F = pQ 2 ^/(vasin\\ (7) 

* E. O. Street, loc. cit.^ p. 352. 
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where p is the density of the water* and v is the kinematic coefficient of 
viscosity appropriate to the circumstances. 

For places on the equator w is zero and equations (3) reduce to 

dz 2 v gx' . dz 2 v dy' 

so that, with the same surface and bottom conditions as before, we have in 
place of (5) 



U Y 2s 



z 2 



Uo V h h 2 ' 
Hence we obtain in place of equation (7) 

3 h ' W 

Now in deep w 7 ater F' is very small compared with pQ 2 >s / (vfl). Therefore 
in estimating the mean rate at which energy is dissipated by the permanent 
currents of all the oceans we may without appreciable error suppose that 
equation (7) is valid everywhere. 

In order to obtain a numerical estimate for the total actual rate of energy 
dissipation due to these non-periodic currents, the surface of the sea has been 
divided into 10-degree " squares" and for each a mean value of Q has been 
used. These values have been deduced from the Admiralty current charts. 
At any point the speed of the current is usually variable ; for the non-periodic 
part, the mean of the two extreme values has been taken. Thus when the 
current is recorded as 20-40 sea miles per day the value 30 has been used. 
The variations in current strengths are greatest in the neighbourhood of the 
equator; owing to the presence of the factor sin X in the formula for the rate 
of energy dissipation the effect of this variation is not important. In any 
10-degree "square" the mean latitude is taken as sufficiently accurate for 
our present purpose. 

The details of the arithmetic are omitted and only the final results for each 
of the three great oceans given. 

For the Atlantic, Pacific, and Indian Oceans the values 

(7, 11, 7) x 10 17 ergs per second, (9) 

have been obtained by substituting in the formula (7), and assuming that 
v = 0*018, the areas over which the summation has been carried out being 

respectively 

(7*3, 13*5, 5-4) x 10 17 sq. cm. 

* In this portion of the work no appreciable error is introduced by regarding p as 
constant, although in all probability the variations of density actually constitute a 
primary cause of the currents. 
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The areas of the surfaces of these oceans are 

(8*2, 16*6, 7*3) x 10 17 sq. cm. 

but there is little information available beyond the fiftieth parallel of latitude, 
or the sixtieth in the case of the North Atlantic. Hence we may say that, 
on the explicit assumption of non-turbulent flow, the presence of viscosity 
causes the non-periodic parts of the ocean currents to produce a dissipation 
of energy at the mean rate 

3 x 10 18 ergs per second. (10) 

It may be noted that this result is roughly 1 erg per second per square 
centimetre of surface. The estimate (10) is to be regarded correct to one 
figure only. 

It is interesting to compare this result with the expression, 22 x 10 18 ergs 
per second, recently obtained by H. Jeffreys* for the mean rate of energy dissi- 
pation in the semi-diurnal tidal motion of all the shallow seas of the world. 
The method there used is due to G-. I. Taylor, f and is based on the assumption 
that the dissipation is caused by direct friction of moving water on the bed 
of the sea, the velocity of slip being taken equal to the surface current 
velocity. The rate of energy dissipation is assumed to be expressible in the 
form| kpv 2 per unit area of bottom, p being the density, v the velocity of slip, 
and k a constant which is taken as equal to that found experimentally for the 
wind blowing over grass land,§ 

The Relations between Currents, Salinities, and Temperatures. 

An estimate of the magnitude of ocean currents can be obtained by simple 
methods from a knowledge of the distributions of the salinity or of the tem- 
perature of the sea water. 

Let s denote the salinity at any point at time t, that is, the mass of salt 

dissolved in unit mass of the water which at that instant surrounds the 

point considered. Making the usual assumption that each small portion of 

water retains its salinity, we have the equation 

.p. 

Wt J 

* 4 Phil. Trans.,' A, vol. 221, p. 239 (1920). 

t 'Phil. Trans., 5 A, vol. 220, p. 1 (1919). 

J Note that this is not of the same form as F or F in equations (7) and (8). 

§ There seems to be no assurance that the agreement between this somewhat arbitrary 
determination of k and another estimate of its value obtained from an extrapolation of 
Bazin's formula for the flow of water in a river may not be largely accidental. The 
various hydraulic formulae for such motion seem to be empirical ; and others, which 
agree equally well with the observations, when extrapolated give results of entirely 
different orders of magnitude, Cf. Lord Eayleigh, { Sci. Papers,' vol. 6, p. 602. 
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where the differentiation follows the path of the water particle. If this 
equation is applied to a particle at the surface of the sea, it can be written in 
the form 

-~ + Qs'cos% + «?^- = 0, (11) 

where Q is the horizontal and %o the vertical velocity of the water, the axis 

of z being measured upwards, s' is the maximum surface gradient of s at 

time t 9 and % is the angle between the directions of Q and s'. 

Similarly, if p is the density and T the temperature of the water at the 

point at time t, the assumption that the conduction of heat is negligible leads 

to the equation 

D (j>T)/Bt = 0, 

while the equation of conservation of mass is 

Bp/Bt = 0. 

From these we deduce 

BT/Bt = 0, 

which as before leads to a result of the form 

~ + QT' cos ^ + w §p = 0, (12) 

Co vZ 

where T' is the maximum surface gradient of T at time t, and ^ is the angle 
between the directions of Q and T'. 

The value of Q could be found from either of the equations (11), (12), if 
the necessary related sets of observations of the salinity or of the tem- 
perature were available. At present, it is possible only to indicate the 
orders of magnitude of the several terms in the equations. It is probable 
that each of the three variables (current, salinity, temperature), is the sum of 
a number of periodic functions of the time, together with a part independent 
of t. The equations should therefore be applied to each of these component 
parts ; but, in view of the scanty nature of the data, this is quite 
impracticable. 

The equations have been applied numerically to the North Atlantic, 
taken as a whole, rather than to particular points in it, the aim, as stated 
above, being to obtain orders of magnitude only, and not definite numerical 
results. A very large number of observations of salinity and of temperature, 
spread over a period of two years, have been collected by H. N. Dickson* ; 
and, from these, he has constructed charts showing the isohalines and the 
isotherms for each of the twenty-four months of the period covered. Their 
most striking feature is the difference between charts relating to the same 

* 'Phil. Trans./ A, vol. 196, p. 66 (1901). 
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month of two consecutive years. It is probable that this variation can be 
attributed to meteorological causes; for it is clear that such events as the 
breaking up of the polar ice will have a profound effect on the chemical and 
physical conditions in the Ocean. 

Adopting the usual practice of measuring the salinity of the water by the 
number of parts by weight of salts dissolved in a thousand parts of water, 
we find as rough mean values* : — 

s' ( = ) 0*2 per 100 sea miles, 

T' ( = ) 2° 0. per 100 sea miles, 
or using C.G.S. units, 

5 '( = )3xl0- 10 s, T'( = )2xl0- 7 . (13) 

In order to obtain the time-gradients of the salinity from the charts, three 
representative points in latitude 50° N"., with longitude 20°, 30°, 40° W., were 
taken, and for each the salinity-time curve was constructed by interpolation. 
The graphs naturally show some irregularities ; but it is evident that the 
predominant parts have a period of six months, their amplitudes being about 
0*1 for the first two stations, and 0*4 for the third. Taking 0*2 as the mean 
amplitude, we find for this term 

ds/dt( = )2xl0~ 9 s. (14) 

The effect of the very small daily variation in salinity is small compared 
with this, notwithstanding the more rapid rate of variation. 

In a similar manner, the time-gradients of the temperature can be 
obtained. As could have been foreseen, the predominating term has a period 
of a year, its amplitude being about 3° C. in each case. This is approximately 
half the mean of the range of surface temperature of all the oceans. There 
is, in addition, a daily component, whose amplitude is about 0*3° C. This is, 
therefore, the more important part, and we may suppose 

8T/8*( = )10- 5 . (15) 

The vertical gradients of the salinity and the temperature are less easy to 
obtain. Thoulett has calculated that the mean value of the salinity of the 
North Atlantic for depths less than 50 m. is 35*82, and for depths from 50 m. 
to 100 m. is 35*99. Using these figures, we find ds/dz( = )Sx 10*" 5 . Other 
observations confirm the accuracy of this estimate. Thus we may take 

ds/dz( = )10-«s, (16) 

* The symbol ( = ) is used to denote equality of orders of magnitude. 

t * Resultats des Campagnes scientif. du Prince de Monaco,' vol. 29, p. 95. 
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while -from a consideration of the vertical variation of temperature at many 

stations we infer 

3T/a»( = )10- 4 . (17) 

The most important part of the vertical velocity of the water particles at 
the surface of the ocean at points far removed from land may be supposed 
due to a semi-diurnal wave of amplitude of the order 30 cm. Hence 

w( = )5xl0- 3 . (18) 

From equations (13), (14), (16), (18) we see that the terms of (11) are of 

orders 

2 x 10- 9 s, 8Q x lO" 1 ^, 5 x 10~ 9 s, 

which are in the ratio 

20, 3Q, 50. 

Thus Q is of the order 20 cm. per sec. or, say, half a knot. 

In a similar manner we see from (13), (15), (17), (18) that the terms of (12) 

&re of orders 

10~ 5 , 2QX10" 7 , 5xl0" 7 , 

which are in the ratio 

100, 2Q, 5. 

Thus Q is of the order 50 cm. per sec. or one knot. 

The mean current velocity of the surface of the Atlantic Ocean is rather 
less than a knot. The agreement of these two rough estimates with the 
observed result is satisfactory, and seems to confirm the mutual consistency 
.of the data employed. 



